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An effort was conducted to numerically determine, using the Nonlinear 
Schrodinger Split-Step Fourier method, if using cross phase modulation could 
cause temporal soliton pulse propagation in a silicon slow-light photonic crys- 
tal waveguide shorter than a millimeter. The simulations demonstrated that 
due to the higher powers and shorter scales of photonic crystals, two-photon 
absorption would cause an optical soliton pulse to be extremely dissipative. 
The model demonstrated, however, that by utilizing cross-phase modulation, 
it is possible to sustain a compressed soliton pulse within a silicon photonic 
crystal waveguide subjected to two-photon absorption over longer relative 
distances. © 2013 Optical Society of America 

OCIS codes: 350.4238, 190.4720. 
1. Introduction 

Optical solitons are undistorted standing waves that result as part of the interplay between 
the nonlinear Kerr effect and anomalous dispersion [1,2]. One method of studying solitons 
numerically is to use the Split-Step Fourier Nonlinear Schrodinger Equation (NLSE) [1-4]. 
The NLSE is one of the most effective methods for numerically analyzing soliton pulse 
propagation in a nonlinear dispersive waveguide. 

Photonic crystal waveguides (PhCWG) have extraordinary nonlinear effects [5], where 
measurements of the dispersion have shown an increase of five orders of magnitude in the 



Group Velocity Dispersion (GVD) compared to the material dispersion [6,7]. PhCWG, with 
their vastly increased nonlinear effects, are an ideal medium for reducing the nonlinear length 
scales from meters and kilometers typical of fibers, to sub-millimeter PhCWG. A soliton is 
formed by the interplay of anomalous GVD, as well as the 3 rd -order Kerr nonlinear optical 
effect. When the pulse is focused by its own intensity, it is known as Self Phase Modulation 
(SPM); whereas when it is focused by a second pulse, it is referred to as Cross Phase Mod- 
ulation (XPM) [1,2]. In addition to large GVD, PhCWG's can be fabricated out of highly 
3 rd -order-nonlinear dielectrics, dramatically reducing nonlinear length scales. 

The first such experimental measurement of soliton pulse compression in a PhCWG was 
conducted using Gallium Indium Phosphate (GalnP) PhCWG [7], which have the advan- 
tage of not having intensity-dependent two-photon absorption (TPA) and TPA-induced 
free-carrier absorption (FCA) [8]. The applications of this technology are numerous, as 
two-dimensional PhCWG may be fabricated into optical waveguides on a device the size 
of a computer chip, to be utilized for all-optical computing [9]. Optical computing has the 
advantage of significantly greater bandwidth compared to current solid-state devices, with 
reduced power and cooling, which can allow vastly improved computational performance [10]. 

Another highly nonlinear dielectric material is silicon [11, 12], which may also potentially 
be used to fabricate chip-scale optical waveguides. Silicon is advantageous, as it is cheap, 
durable, readily available, and there already exists an infrastructure to mass-produce silicon 
devices. Pulse compression was previously observed experimentally in novel silicon nanowires 
[13,14], but previous analytical and numerical studies have demonstrated that with SPM 
alone, TPA and FCA will prevent a pulse from self-focusing itself sufficiently to generate a 
soliton-like wave [8]. Even if the soliton pulse were to compress, it would simply dissipate 
after propagating relatively short nonlinear lengths. One of the advantages of XPM is that a 
pulse can be focused without needing to be highly energetic and intense [1]; therefore there 
was interest in seeing if a second pump pulse could be used to focus the primary signal pulse, 
to induce the soliton. If XPM can be used to overcome the limitations of TPA and FCA 
and enable chip-scale soliton propagation in silicon, it may very well develop into practical 
all-optical computing. 

2. Nonlinear Schrodinger Equation 

In order to accurately model the optical soliton, two crucial nonlinear characteristics are 
necessary, the GVD and the nonlinear Kerr effect [1]. The GVD is determined from the 
derivatives of the refractive index over the angular frequency. The other component of soliton 
pulse propagation is the nonlinear Kerr effect, which results in a change of refractive index 
proportional to the square of the electric field, or proportional to the optical intensity. When 
there is strong SPM and/or XPM interplaying with strong anomalous GVD, the SPM/XPM 
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will cause the leading edge frequencies to be lowered, and the trailing edge frequencies to be 
raised. At the same time, the anomalous GVD will cause the lowered leading edge frequencies 
to slow down, and the raised trailing edge frequencies to speed up. This results in the pulse 
compressing temporally into a soliton pulse [1,2]. 

The GVD is determined by knowing the derivatives of the refractive index over the angular 
frequency. These terms are found using the following equations [1,2,4]: 
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where n(u) is the refractive index as a function of angular velocity, u (rad/s) is the angular 
frequency, c (m/s) is the speed of light, vq (m/s) is the group velocity, and no is the group 
index. These values can be determined either experimentally [11,15], approximated with the 
Sellmeier equation [1,2], or engineered to specific values with geometric dispersion, as is the 
case of photonic crystals [5,6]. 

The other component of generating soliton pulse propagation and compression is the non- 
linear Kerr effect. The change in refractive index is determined by: 

An = n 2 -(J 1 + 2-/ 2 ), (5) 

1 = (6) 

where n 2 (m 2 /W) is the nonlinear Kerr coefficient, r] is the vacuum impedance (~ 1207T 
Q), n is the refractive index, A (V/m) is the electric field, and l\ and J 2 (W/m 2 ) are, 
respectively, the intensities of the pulse and another coupled pulse. This nonlinear term 
can be determined from the 3 rd -order nonlinear susceptibility, which is also responsible for 
four- wave mixing [2,16]: 

x (3) 
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where eo is the electric permittivity in a vacuum (8.854 -10 12 F/m), and (m 2 /V 2 ) is the 
3 rd -order nonlinear susceptibility coefficient [2]. Throughout this study, a Kerr nonlinearity 
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of 4-10 18 m 2 /W [6,11,17] is used. The nonlinear phase shift is thus [1]: 

A(f) NL = n 2 -k -Lih + 2-h) = n 2 - — -L-{h + 2-/ 2 ), (8) 



where ko (m -1 ) is the optical wave- number, A (m) is the wavelength, and L (m) is the optical 
propagation length. These two nonlinear optical effects can be combined to form the NLSE, 
which in its complete form is represented as: 

dA i n dA 2 i n dA 3 . IAl2A 

where (3 2 (s 2 /m) is the 2 nd -order GVD coefficient, @z (s 3 /m) is the 3 rd -order GVD coefficient, 
A e (m 2 ) is the effective cross-sectional area of the waveguide, 7 (m-W) -1 is the nonlinear 
parameter, a (m _1 ) is the linear loss coefficient, i = \J — 1, and A (V/m) is the pulse electric 
field envelope function. 

The NLSE is also known as the Split-Step Fourier method because each incremental dis- 
tance step involves two fundamental steps. The first step is to account for the GVD. A code 
will convert the pulse information into the spectral domain by using the Fast Fourier Trans- 
forms (FFT) method [1,4,18,19], apply the GVD, and then convert the pulse back from the 
spectral domain to the temporal domain to calculate the Kerr nonlinearity. The equation for 
the GVD is as follows [1,4]: 
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where A(z,T) and A(z,u) are the pulse envelope functions in the temporal and spectral 
domains, T is the normalized time, v g is the group velocity, z is the propagation distance, 
and lj is the normalized angular frequency. The second step is to take the pulse, now in the 
temporal domain, and apply the effects of SPM/XPM. The equation for the nonlinear phase 
shift is [1,4]: 
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After completely these two steps, one can apply the effects of linear and nonlinear loss for 
the specified distance increment. By following this algorithm, soliton pulse propagation can 
be simulated and studied. 



3. Two Photon Absorption 

One big challenge to soliton pulse propgation in silicon is the nonlinear effect of multi- 
photon absorption [8]. Dielectric materials have an electronic band gap, where photons with 
energy greater than the electronic band gap get absorbed by the material [20] . The electronic 
band gap for silicon is 1.1 electron-volts (eV) [7,8]. As it is desired to use light pulses at a 
wavelength of 1550 nm, this corresponds to a photon energy of 0.8 eV, and thus silicon is 
transparent at this wavelength. Multi-photon absorption, however, is the result of multiple 
photons interacting, and their combined energy exceeding the electronic band gap. For ex- 
ample, at 1550 nm, two photons will result in a total energy of 0.8x2=1.6 eV, which exceeds 
the 1.1 eV electronic band gap of silicon. Therefore, silicon is subject to TPA at 1550 nm. 

In order to calculate the impacts of TPA on silicon, one must find the coefficient of TPA 
at the given wavelength by the following equations [21]: 
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where EP is equal to 21, h is Plank's constant over 2tt, v is the frequency, no is the refractive 
index, and E g is equal to the 1.1 eV electronic band-gap of silicon [7,8,12]. The constant 
if was calculated based on the experimentally measured value of 4.4-10~ 12 (m/W) at 1550 
nm [22]. With the value of /3 T pa, the TPA can be simulated with [12,23]: 

I(z,t) = I(0,t)-exp[-PT PA -I{0,t)-z], (18) 



Once the two-photon absorption coefficient is determined, one can determine the impact 
from the free carriers generated by the TPA. For ultrashort pulses, where the pulse is sub- 
stantially shorter than the free-carrier lifetime, the FCA tends to be substantially less than 
the effects of TPA [12]. The model calculates the FCA at each distance step separately, 
assuming that there were no free-carriers there at the beginning. This is an approximation, 
as carriers will in fact drift over to the different distance steps, but this is difficult to ac- 
curately predict. For this reason, the assumption that carriers stay localized will be used. 
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After the SPM/XPM and dispersion was applied to the pulse, the model predicts the rate 
of free-carrier generation from the two-photon absorption [12,22-24], 



dN(t) 



= G - 



N(t) 



(19) 



dt 
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where (3tpa (m/W) is the coefficient of two photon absorption, I (W/m 2 ) is the optical 
intensity (either from the signal or coupled pulse), v is the frequency, H is Planks Constant 
over 2n, G is the rate of free-carrier generation, r TPA (s) is the free carrier lifetime, and N(t) is 
the free-carrier density. Once the rate of free-carrier generation over time has been determined 
at the specific distance step, the next step is to calculate the density of free-carriers. While 
the simulation goes through the process of generating carriers over time, it also reduces 
carriers based on the specified free-carrier lifetime. After the carrier density is determined as 
a function of time (at the given distance increment), the appropriate free-carrier absorption 
is applied to the pulse in the temporal domain. The equation is as follows [12,23]: 



By implementing these equations within the NLSE simulation, an accurate analysis of soliton 
pulse propagation in a silicon waveguide subjected to TPA and FCA can be achieved. 

4. Simulations of Cross Phase Modulation 

The NLSE model that will be used throughout this study was first set to a single pulse with 
SPM. First, the model propagated a 10 pJ hyperbolic secant squared pulse down a 10 cm 
silicon nanowire; the results of this simulation are shown in Figure 1-a. A compression ratio 
factor of 7 was clearly observed. The model then scaled the pulse energy and waveguide 
length to match the soliton number [7] and nonlinear length [1] as the nanowire study, using 
the following relations: 



where L (m) is the waveguide length, E (J) is the input pulse energy, /3 2 (s 2 /m) is the 2 nd - 
order GVD coefficient, A e (m 2 ) is the waveguide effective cross-sectional area, and ric and 
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n are the group and refractive indices, respectively. Based on these relations, it is clear that 
by using a PhCWG, the soliton length scales can become reduced dramatically. The caveat 
is that it is necessary to use higher pulse-energies, which in turn have greater TPA and 
FCA, which attenuate the soliton, and disturb the pulse, as observed in Figure lb. While 
compression has been observed at the early onset of soliton propagation within the PhCWG, 
these compressed pulses simply cannot sustain themselves in the presence of intense TPA 
and FCA. 

The NLSE model was modified to study both XPM as well as the effects of cross-TPA 
(XTPA) and the resulting FCA [12]. In this XPM simulation the process is completed twice, 
first for the primary pulse, and second for the secondary pulse. It uses a modified split-step, 
since there are two coupled pulses, to take into account differences in the group velocity 
which result from the different wavelengths between the pump and the signal. 

Part of the modified aspect is a loop for the nonlinear effects that require considering 
both pulses, the two include XPM, XTPA, and FCA induced from XTPA. The split-step 
calculations are used in a normalized time domain [1], however, with different wavelengths 
and thus different group velocities, this normalized time is not always applicable. In fact, 
one concern with XPM is the case of walk-off where a faster traveling pulse fully clears away 
from the slower pulse within the waveguide [25] . A time shift factor is determined, based on 
the difference in group velocities, defined as 5T, where 

ST = (n'c-n G )iJ^), (24) 

where Sz is the distance increment, 5t is the time step used in the simulation, and Co is the 
speed of light in a vacuum. The value of ST is rounded down, as it is an integer number of 
time steps to shift, and added to calculate the effects of pulse 2 onto pulse 1 (and subtracted 
for pulse 1 onto pulse 2). 

5. Soliton Propagation with Cross Phase Modulation 

This model was run using the parameters experimentally observed in silicon PhCWG [6], 
with a negligibly low pulse energy (to avoid SPM) of 10 fJ pulse for the signal pulse 
near the highly-dispersive photonic band-gap wavelength at 1544 nm, coupled temporally 
with an intense 105 pJ pump pulse at 1536 nm. In this particular example, the GVD at 
the signal wavelength was a highly dispersive anomalous -4 ps 2 /mm, whereas the pump 
wavelength had a normal GVD coefficient of 100 ps 2 /m. The signal was also set to a 
third-order GVD coefficient of 0.1 ps 3 /mm, a value typical of these engineered photonic 
crystal waveguides [7]. While higher-order dispersion is known to affect soliton propagation, 
especially at wavelengths very close to the highly dispersive photonic crystal band-gap [26], 
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there was negligible difference in outcomes for the practical lengths modeled in this study. 
Finally, a linear loss of 2.4 dB/mm was added to the simulation. 



By the phenomena of XPM, temporal compression was observed; the signal pulse under- 
went compression by a factor of 5 after propagating through 120 /im of the PhCWG. This 
result is displayed both as a function of normalized time and propagation distance in Figure 
2-a; as well as a comparison of the normalized temporal intensity function in Figure 2-b. 
Finally, the pump pulse did not experience any compression, as it was subjected to much 
less GVD; this is demonstrated in Figure 2-c. As a result of this, the simulation clearly 
demonstrates that the compression is the result of XPM, as no compression was observed 
when the same 10 f J signal pulse propagates through the waveguide without the pump pulse. 

When running the simulation, a numerical scale was developed to represent XPM-induced 
solitons. Much like a SPM-soliton [7, 8] , the XPM-soliton number N is a ratio of the funda- 
mental soliton energy: 

» = /f • < 25 > 

E, = 2 M, (26) 



where 02 (s 2 /m) is the 2 nd -order GVD coefficient, 7 (W-m) -1 is the nonlinear Kerr coefficient, 
r (s) is the temporal duration of the signal pulse in seconds, and E p (J) is the pump pulse 
energy. The pulse compression ratio requires higher XPM soliton numbers than with SPM; 
this is expected due to the effect of pulse walk-off, which becomes more profound at longer 
waveguide lengths. At the practical length scale of 0.1 nonlinear lengths (similar length 
studied with SPM), which equated to 0.4 mm (a typical photonic crystal waveguide [6-8, 
11]), a pulse compression of nearly 3 was observed from an XPM-soliton number of 15; the 
results of these studies can be seen in Figure 3-a. The Time Bandwidth Product (TBP) was 
tabulated and is represented by Figure 3-b. Under the practical conditions mentioned earlier, 
the TBP remained relatively unchanged; evidence of XPM-induced soliton pulse propagation. 
At longer length scales, the pulse eventually underwent broadening and the TBP increased; 
this can be attributed to TPB attenuating the pump pulse, as well as pulse walk-off. The 
nonlinear loss (in dB) was also observed in figure 3-c. In all cases, the longer the nonlinear 
propagation length, the more loss was observed; this is to be expected. It was interesting 
to note, that in shorter propagation lengths, higher XPM-Soliton numbers had an ability to 
actually increase the peak energy as the signal pulse was compressed by XPM. In all cases 
studied, the optical pulses were all by-definition dissipative solitons; it was demonstrated 
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that XPM can sustain the pulse over much longer propagation lengths scales than SPM 
alone. 

6. Practical Applications of Cross Phase Modulation 

The previous study was on the assumption that both the pump and signal wavelengths 
were subjected to TPA and FCA; this may be practically realized by utilizing a spectrally 
broad pulse and narrow bandpass filters [25]. While this may be simple to implement, it has 
the limitation of significant power attenuation from the filters. Another potential option to 
generate optical pulses of different wavelengths coupled temporally and spatially is to use an 
Optical Parametric Oscillator (OPO), where a broad range of wavelengths may be generated 
by adjusting the length of the resonating cavity [27,28]. In an OPO, signal and idler pulses 
are generated with a combined wavelength energy equal to the pump wavelength energy: 

111 

A = A +A - • (27) 

^pump ^signal ^idler 



By utilizing this technology, one may realistically generate pulses coupled temporally and 
spectrally, where the pump pulse (not to be confused with the pump pulse for the OPO) has 
a wavelength not subjected to TPA within silicon. In this next study, the authors investigate 
pulse compression resulting from a low-powered signal pulse at 1550 nm; and a pump pulse at 
a wavelength of 2300 nm, which is not subjected to TPA due to silicon's electronic band gap. 
This can be achieved with a properly aligned OPO pumped 926 nm, a realistic wavelength 
that can be achieved in a titanium sapphire laser [29]. 

The study was conducted with the same dispersion and energy parameters as before, but 
without the pump pulse subjected to TPA, and a compression ratio over 5 was observed 
at 90 fim of propagation. This is shorter than the 120 /xm of propagation needed for the 
study in Figure 2; and the results of this study can be observed both in Figure 4-a and 
Figure 4-b. It was interesting to observe that the generation of free carriers, as shown in 
Figure 4-c, actually increased as the signal pulse, which is subjected to TPA, compressed 
and increased in temporal intensity due to the XPM from the pump pulse. In addition, it 
was observed that the compression remained consistent for even the longer propagation 
studies. This study was repeated for a host of different XPM Soliton Numbers (Figure 5-a) 
and nonlinear length ratios (Figure 5-b), and the nonlinear losses were observed (Figure 
5-c). These nonlinear losses were very close to the previous non-OPO study (Figure 3-c). 
The data clearly shows that by utilizing a substantially longer wavelength pulse, practical 
soliton pulse propagation in silicon can be both achieved and sustained for low-energy 
pulses at wavelengths subjected to intensity-dependent TPA. 
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7. Conclusion 

The simulations conducted have demonstrated the physical ability of soliton pulse com- 
pression and propagation via XPM in a silicon PhCWG with TPA and FCA present 
at both the pump and signal wavelengths. Use of a pump pulse and XPM to focus the 
pulse enables the soliton phenomena to occur on pulses with much lower intensities and 
thus much lower TPA. One challenge that was successfully overcome was to develop a 
NLSE Split-Step Fourier algorithm that would take into account the differences in group 
velocity and pulse walk-off that will occur between two pulses of different wavelengths. The 
simulation demonstrated that through the use of XPM, temporal soliton pulse propagation 
can be realized over practical length scales with low-intensity ultrashort pulses in a silicon 
PhCWG with TPA. This phenomenon might be practically used to sustain a data pulse 
within a silicon waveguide, a capability with numerous potential technological applications. 
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Fig. 1. SPM-induced soliton propagation, in a (a) silicon nanowire and a (b) 
silicon PhCWG. 



Primary Pulse - 1 544 nm Waveguide Output - XPM Waveguide Output - No XPM 




Fig. 2. XPM induced compression: (a) is the signal pulse compression as func- 
tion of waveguide propagation distance and time; (b) is the waveguide output 
of weak signal pulse under XPM; and (c) is the waveguide output of weak signal 
pulse without XPM, thus proving the compression is the result of XPM. 
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Fig. 3. Results of numerical XPM studies: (a) the signal pulse compression 
factor, (b) the output signal pulse TBP, and (c) the nonlinear peak intensity 
loss (dB). 




Fig. 4. Results of numerical XPM studies, with a low-energy C-band pulse cou- 
pled with a high-energy 2300 nm pump pulse: (a) signal pulse compression as 
function of waveguide propagation distance and time; (b) numerical predicted 
free-carrier density as a function of waveguide propagation distance and time; 
and (c) waveguide output of weak signal pulse under XPM. 
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Fig. 5. Results of numerical XPM studies, with a low-energy C-band pulse cou- 
pled with a high-energy 2300 nm pump pulse: (a) the signal pulse compression 
factor, (b) the output signal pulse TBP, and (c) the nonlinear peak intensity 
loss (dB). 
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